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According to a classical result due to Hudson, the Wigner function of a pure, continuous variable 
quantum state is non-negative if and only if the state is Gaussian. We have proven an analogous 
statement for fìnite-dimensional quantum systems. In this context, the role of Gaussian states is 
taken on by stabiUzer states. The general resuhs have been pubUshed in [D. Gross, J. Math. Phys. 
47, 122107 (2006)]. For the case of systems of odd prime dimension, a greatly simplified proof can 
be employed which stili exhibits the main ideas. The present paper gives a self-contained account 
of these methods. 



I. INTRODUCTION 

The Wigner distribution establishes a correspondence 
between quantum mechanical states and real pseudo- 
probability distributions on phase space. 'Pseudo' refers 
to the fact that, while the Wigner function resembles 
many of the propcrties of probability distributions, it 
can take on negative values. It is therefore of interest 
to charactcrizc thosc quantum states that are classical 
in the sense of giving rise to non-negative phase space 
distributions. 

For the case of pure states described by vectors in 
H = L^(R), the resolution of this problem was given 
by Hudson in Ref. and later extended to multiple- 
particles by Soto and Claverie (Ref. Q)- 

Theorem 1. f Hudson, Soto, Claverie) Let i' e ^^(R") 
be a state vector. The Wigner function of ip is non- 
negative if and only if ìp is a Gaussian state. 

By definition, a vector is Gaussian if and only if it is 
of the form 

where x G R" and 6 is a symmetric matrix with entries 
in C !M]- 

It is our objective to prove that the situation for dis- 
crete quantum systems is very similar, at least when the 
dimension of the Hilbert space is odd. The following 
Theorem states the main result. 

Theorem 2. (Discrete Hudson's Theorem) Let d he odd 
and ip G L'^{7L'^) he a state vector. The Wigner function 
of ìp is non-negative if and only if tp is a stabilizer state. 

Given that ìp{q) ^ for ali q, a vector ìp is a stabilizer 
state if and only if it is of the form 

ìpiq) = e^'f-J^^+^-j), 

where x G "Z^ ^'^'^ 9 is a symmetric matrix with entries 
in 'Ed- 
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In the previous theorem, := {0, . . . , — 1} denotes 
the set of integers modulo d. 

It turns out that, although the formulation of Hudson's 
result carries naturally over to finite dimensionai systems, 
the respective proofs are radically different. The originai 
argumcnt relies crucially on function theory, which is of 
course not available in the setting that this paper ad- 
dresses. 

Recently, Galvao et. al. took a first step info the di- 
rection of classifying the quantum states with positive 
Wigner function (see Ref. [3[). To explain the relation- 
ship of thcir rcsults to the present paper, we have to 
commcnt shortly on two different approaches to defin- 
ing discrete Wigner functions. On the one hand, it 
has long been realized that Wigner's definition carries 
over naturally to discrete odd-dimensional systems (Ref. 
3, i, S S S S [S ini, m, [H, H and Sectlon |TTJ) . This 
approach is the one used in the present paper. On the 
othcr hand, Gibbons, Hoffmann, and Wootters listed a 
set of axioms that candidate definitions have to fulfiU in 
order to resemble the properties of the well-known con- 
tinuous case (Ref. [l^l)- Let us cali functions that fall 
info this class generalized Wigner functions. The char- 
acterization does not specify a unique solution: for a ri- 
dimensionai Hilbert space, there exist d*^"^ distinct gen- 
eralized Wigner functions. The construction of Gibbons 
et. al. has been described only for the case where d is 
the power of a prime. 

If the dimension of the Hilbert space is of the form c?" , 
a second ambiguity arises. We are free to conceive such 
a space either as being associated to a system of n con- 
stituents, each of dimension d, or to a single one of dimen- 
sion d". While the Wigner function is the same for both 
cases, the set of stabilizer states is not (see Refs. (sl.!!^). 
Indeed, the 'single-particlc' stabilizer states turn out to 
be a proper subset of the 'multiplc-particle' oncs. As 
a striking examplc, the generalized Bell and GHZ states 
nor are not stabilizer states on single or d^-dimensional 
systems. 

In Ref. 0] it was proved that a state of a single-particle 
system of prime-power dimension is a stabilizer state if 
and only if ali its generalized Wigner functions are non- 
negative. The authors aim to establish necessary require- 
ments for quantum computational speedup. Indeed, if 
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the Wigner function of a quantum computer is positive 
at ali times, then it operates only with stabilizer states 
and hence offers no advantage over classical computerà, 
by the Gottesman-Knill Theorem (Ref. [iBl)- 

For the case of pure states, our results imply the ones 
of Ref. 3 and exceed them in two ways. Firstly, we show 
that it suffices to check positivity for a single definition of 
the Wigner function, as opposed to d'^~^ ones. Secondly, 
our statements hold for multiple-particle systems, which 
constitute the proper setting for both quantum compu- 
tation and the Gottesman-Knill Theorem. On the other 
hand, Ref. 0] makes assertions about mixed states and 
qubit systems, which are not covered by our findings. 

Our general results have been published in Ref. 
However, the proof is rather involved. Many technicali- 
ties arise due to the fact that for non-prime d, arithmetic 
modulo d lacks the desirable properties of finite fields. 
Our aim in writing Ref. Q was to achieve the broad- 
est possible generality in spite of these difficulties. The 
downside of this approach is that core ideas of the argu- 
ment are obscured by technical issues. The present paper 
employs a different method of proof, which is available 
only for systems of odd prime dimension. For this spe- 
cial case, the main result of Ref. 0] can be obtained using 
only a fraction of the space. It is our hope that this paper 
makes the ideas accessible to a wider audience. 

The next scction summarizcs furthcr findings con- 
tained in Ref. We go on to recali the definition and 
properties of discrete Wigner functions in Section Hill 
Section [Tv] is devoted to a complete proof of the easiest 
special case of Theorem [2l that being given by a single 
particlc on a Hilbert space of prime dimension. 



II. FURTHER RESULTS AND IMPLICATIONS 

It is naturai to ask how Hudson's results generalize to 
mixed states. Certainly, mixtures of Gaussian states are 
positive on phase space and Narcowich in Ref. con- 
jectured that ali such quantum states are convex combi- 
nations of Gaussian ones. Bròckcr and Werner refuted 
the conjecture by giving a counter-example (Ref. [l3l)- 
We show in Ref. ^ that the situation is similar in the 
finite setting. 

Further, we show how to lift the ambiguity in the ax- 
iomatic characterization of Wigner functions by requir- 
ing Clifford covariance, note that a unitary operator pre- 
serves positivity if and only it is a Clifford operation, 
discuss the relation of various ways to introduce Wigner 
functions and stabilizer states in dimensions of the forni 
d = p", and give an explicit account on the connection 
between stabilizer states and Gaussian states. 



III. WIGNER FUNCTIONS 

This section provides a very superficial introduction 
to discrete Wigner functions. Wc allow ourselves to refer 



the reader to Refs. ^ for further details. In what follows 
d denotes an odd prime. AH integer arithmetic in this 
paper is implicitly assumed to be modd. The symbol 
2^^ = (o?-|-l)/2 is the multiplicative inverse modulo d. AH 
state vectors are elements of the Hilbert space Ti spanned 
by {|0), . . . , |d — 1)}. Lastly, uj ~ e^* is a dth root of 
unity. 

The relations 

x{q)\k) ^\k + q), zip)\k) = UjP''\k) 

define the shift and boost operators respectively. The 
most centrai element in the theory are the Weyl operators 
(in quantum Information also known as the generalized 
Pauli operators) given by 

w(j),q)~uj^^ ^'^z(p)x{q). 

The characteristic function of an operator p is given by 
the expansion coefficients of p in terms of the Weyl op- 
erators 

We define the Wigner function to be the symplectic 
Fourier transform of the characteristic function: 

Wp^P'l) = ^ E ^^^"'"S,(?,x). 

It is a tedious yet straight-forward computation to show 
that the Wigner function of a pure state is given by 

Wi,{p,q) Wi^,)(^^,i{p,q) 

= ^ E ^"^^^(9 + 2-^0^5(9 -2-ie). 

If 5 is a 2 X 2-matrix with elements in and determi- 
nant 1. then there exists a unitary operation fJ.{S) (the 
Weil [l^l [m or metaplectic representation of S) such 
that 

p{S)w{p,q)p{S)'^ ^w{S{p,q)). 

The Wigner function is covariant in the scuse that, if 
p' = p{S)ppiSy, then 

Wp,ip,q) = Wp{SÌp,q)). (1) 

Similarly, the Weyl operators induce translations of 
the Wigner function. Lctting p' = w{p\q') pw{p\q')\ 
it holds that 

Wp.{p,q)=Wp{p + p\q + q'). (2) 

The Clifford group is the set of unitary matrices that 
send Weyl operators to Weyl operators under conjugation 
[23 |. Every Clifford mapping is of the form w{p,q)p{S) 
and hence preserves positivity of the Wigner function. 

Finally, stabilizer states are the images of the com- 
putational basis states under the action of the Clifford 
group. 
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IV. MAIN THEOREM - SINGLE PARTICLES 
IN PRIME DIMENSIONS 

Define the self correlation function 
and note that the Wigner function obeys 

X 

Recali that the Fourier transform / of a function / : 
^ C is defined to be f{x) = l/dEg^"'^''/('?)- 
Therefore, for a fixcd qo, W{p,qo) is the Fourier trans- 
form of K{qQ, x). Hence W is non- negative if and only if 
the d functions K{qQ, ■) have non-negative Fourier trans- 
forms. 

In harmonic analysis, the set of functions with non- 
negative Fourier transforms is characterized by a well- 
known theorem due to Bochner. We state an elementary 
version of Bochner's Theorem, along with a variation for 
subsequent use. 

Theorem 3. (Variations of Bochner's Theorem) Con- 
sider a function / : ^ C • It holds that 

1. The Fourier transform of f is non-negative if and 
only if the matrix 

A\ = f{x - q) 

is positive semi-definite. 

2. The Fourier transform of f has Constant modulus 
(i.e. \f{x)\ = constj if and only if f is orthogonal 
to its translations: 

(/,Ì(g)/)=^/»/(x-Q) = 0, 

X 

for ali non-zero q ^TLd- 
Proof. The matrix A is circulant. It is wcU-known that 
circulant matrices are normal (hence diagonalizable) with 
eigenvalues given by the Fourier transform of the first row 
(up to a positive normalization Constant). The first claim 
is now immediate. 

By the same argument, A is proportional to a unitary 
matrix if and only if |/(g)| is Constant. But a matrix is 
unitary if and only if its rows form an ortho-normal set 
of vectors. □ 

The next three lemmas harvest some consequences of 
Bochner's Theorem to gain Information on the pointwise 
modulus \4'{<ì)\ of the vector. 

Lemma 4. (Modulus Inequality) Let be a state vector 
with positive Wigner function. 
It holds that 

mq)f > mq - x)\\i:{q + x)\ 

for ali q,x G Z^. 



Proof. Fix a g G Z^. As W.^ is non-negative, so is the 
Fourier transform of K^{q,x) with respect to x. Boch- 
ner's Theorem implies that A'^y ~ K{x — j/, q) is positive 
semi-definite {psd) which in turn implies that ali princi- 
pal sub-matrices are psd. In particular the determinant 
of the 2x2 principal sub-matrix 

( K^{q,0) K^{q,2x)\ 
\K^{q,-2x) K^iq,0) ) 

( _ \^{q)? ìP{q + x)i){q-x) \ 

\il;{q + x)il;{q-x) ) 

must be non-negative. But this means 

|^(<Z)|4- |Vi(g + x)^((7-x)|2 >o, 

which proves the theorem. □ 

We will cali the set of points where a state-vector is 
non-zero its support. 

Lemma 5. (Support Lemma) Let ip be a .state vector 
with positive Wigner function. 

If Tp is supported on two points, then it has maximal 
support. 

Proof. Denote by S* = supp ip the support of ìp. S has 
the property to contain the midpoint of any two of its 
elements. Indeed, if a, ò G S, then setting q ~ 2~^{a -\- b) 
and X = 2~^{a — b) in the Modulus Inequality shows that 

|V(2-i(a + fe))| > |^(a)||V(6)| >0, 

hence 2-'^{a + b) e S. 

Assume there exist two points a, ò G S. Requiring 
a = is no loss of generality, for else we substitute ìp 
by ip' = w{0, ~a)ip. By Eq. ìp' has positive Wigner 
function if and only if ìp has. 

We claim that 

2~'l3b&S (3) 

for ali l and f3 < 2K The proof is by induction on 
Suppose Eq. ^ holds for some l. If /3 < 2'+^ is even, 
then 2-'-i/3ò = 2-'(/3/2)fe G S. Else, 

2-'-i/3ò = 2-i(2-'^fe + 2-'^fe) eS, 

which proves the claim. 

Now,by Fermat's Little Theorem 2'^~^ = 1 mod d and 
hence, setting / = d — 1 in Eq. Q, we conclude that 
Pb E S for ali /3 < d — 1 . But every point in "Ed is of that 
form. □ 

Lemma 6. (Constant Modulus) Let ìp be a state vector 
with positive Wigner function and maximal support. 
Then \ìp((l)\ = const. 
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Proof. Pick two points x,q € 'Ed and suppose > 

Letting z ~ x — q, the assumption reads |'0('?)l > 
2;)|. Lemma m centered at g + z gives 

|^(q + z)|2 > |V;(g)||V(g + 2z)| 

> |7/.(9 + z)||i/;(g + 2z)|, 

therefore \tfj{q + z)\ > \t/j{q + 2z)\. By inducting on this 
scheme, we arrive at 

H{q)\ > Hiq + z)\ > \^{q + 2z)\ > ••• 

and hence \'>p{q)\ > \ìp{q + dz)\ ~ \'ìp{q)\, which is a con- 
tradiction. 

Thus \ip{q)\ < Swapping the roles of x and q 

proves that equality must hold. □ 

Theorem 7. (Main Theorem - Special Case) Let d he 
prime and ip S L'^{7Ld) he a state vector with positive 
Wigner function. Then tp is a stahilizer state. 

Proof. By the Support Lemma, tp is eithcr a position 
eigenstate or else it has maximal support. In the formcr 
case, tp is manifestly a stabilizer state, so we need only 
treat the latter. Let U he a Clifford operation. Since U 
preserves positivity, the Support Lemma applies io Uip. 
Suppose U is such that suppC/-0 contains just a single 
point. Then Uip belongs to the computational basis and 
hence, by definition, ip ìs & stabilizer state. 

Therefore, we are left to treat those state vcctors whose 
image under any Clifford operation has maximal support. 
The proof is concludcd by showing that such statcs do 
not exist. 

For assume there is such a vector tp. As ip has point- 
wise Constant modulus, so does K.^. Employing Theorem 
[31 we find that, for every fixed qo, W{p,qo) is orthogo- 
nal to its own translations. But since W is non-negative, 
it foUows that W{p, qo) can bc non-zero on at most one 
point. A Wigner function that is concentrated at a sin- 
gle point can not rcprcscnt a physical state [2^. There 



must hence exist at least two points a, h in the support 
of W (note that we are now considering the support of 
Wigner functions and no longer the support of state vec- 
tors) . Making once more use of the fact that translations 
are implemented by Clifford operations, assume a = 0. 
There exists a unit-determinant matrix S that sends b to 
a vector of the forni Sb ~ {0,qo)^. But then there are 
two points in the support of W^(s)^(p, qo), contradicting 
our earlier derivation. □ 



V. SUMMARY 

We bave proved a 'classicality result' for discrete Wig- 
ner functions: those state vectors which give rise to a 
classical probability distribution in phase space belong 
to the set of stabilizer states. These, in turn, allow for 
an efficient classical dcscription. Comparing the proof of 
the special case treated bere to the involved argument 
employed in Ref. @, it becomes apparent how much the 
geometrical properties of integer residues modulo prime 
numbcrs simplify the structure. 
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